On the physical origin of fatigue phenomena 
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The fractional power-law material behavior has been investigated within the framework of a 
modified mean field theory, in which high-temperature structure precursors in a crystalline or poly- 
crystalline material are treated as a partially ordered liquid over a wide range of temperature. By 
doing so, the relaxation behavior of the material can be studied via a modified Landau-Khalatnikov 
equation. It then becomes clear that, as the special cases of the fractional power-law behavior, 
seemingly quite different fatigue phenomena and physical deterioration are governed by a funda- 
mental physical phenomenon, i.e., the competition between ordered structures and partially ordered 
high-temperature structure precursors during a series of atomic relaxation processes. 

PACS numbers: 61.43.-j, 62.20.me, 64.60.Cn, 64.60.F- 



I. INTRODUCTION 

Most material properties are often characterized by a 
linear relationship between the applied field and the ma- 
terial response. For example, the stress-strain relation- 
ship of a test specimen within the elastic range can be 
written as a = Ee 1 here a and e represent the applied 
stress and the induced strain, respectively; E is defined 
as the elastic modulus. In many cases, however, such 
a linear relationship could deteriorate into a fractional 
power-law relationship. For instance, the stress-strain 
relationship of the test specimen undergoing plastic de- 
formation might be a = K s ht n (Hollomon's equation [l[), 
here n is defined as the strain hardening exponent and 
K s h is usually called the strength index; in practice, for 
most metallic materials, n generally varies between 0.2 
and 0.5 Q. Such examples could be multiplied indefi- 
nitely. Generally speaking, the fractional power-law be- 
havior is very often connected to physical deterioration 
in solid materials. In one specific field, material fatigue, 
this trend is particularly evident. Over the past two hun- 
dred years, material fatigue has been extensively studied 
and many fatigue laws have been developed. Those laws, 
such as the Coffin-Manson law 0, 0] , the Basquin equa- 
tion [f|, the Wohler equation or the Wohler curve 
etc., are empirical relationships and their validity and 
reliability have been tested and verified by numerous ex- 
perimental data. But their microscopic physical origins 
have not been clarified. It is interesting to notice that 
the Coffin-Manson exponent has shown remarkable uni- 
versality in single-phase metallic materials 0] , which has 
raised a question if there is a fundamental physical phe- 
nomenon underlying apparently dissimilar fatigue phe- 
nomena in different materials. Surprisingly, the above- 
mentioned fatigue laws are also fractional power-law rela- 
tionships, which provides a perfect opportunity to check 
whether there exists a universal phenomenon that gov- 
erns the fractional power-law material behavior. If it does 
exist, then, at least, those fatigue laws should have the 
same origin. In this paper, we will try to give our answer 
to this question. 



II. STATISTICAL INTERPRETATION OF 
DISORDERED STRUCTURES 



It is well known that physical deterioration and fa- 
tigue phenomena in solid materials always involve a large 
amount of disordered microscopic structures. Therefore, 
a good description of those material behavior should be 
able to provide a statistical interpretation to explain the 
emergence and evolution of disordered structures in the 
test specimen undergoing external loading. It is perhaps 
worth briefly mentioning how statistical methods have 
been used to study the occurrence of fracture, which is 
slightly different from our studies. Some methods, such 
as fractals, scaling laws, percolation networks, etc., have 
been exploited to investigate fatigue and fracture phe- 
nomena d, ■ These advanced statistical methods have 
successfully resulted in some theoretical models that can 
be used to predict fracture and cracking behavior in ma- 
terials; the predicted results by these models are in good 
agreement with experimental data to some extent. For 
instance, the percolation theory was used to describe the 
occurrence of fracture as a critical phenomenon Q and 
such a conceptual model was supported by some exper- 
imental results, like that reported in Ref. The 
problem associated with these statistical methods is that 
disordered structures inside materials are neither purely 
stochastic in the time domain nor randomly distributed 
in the spatial domain over a wide range of temperature. 
Thus, those models may not be able to correctly de- 
scribe how disordered structures alter material behav- 
ior and eventually lead to the occurrence of fracture and 
crack. In our opinion, any statistical method that can 
correctly describe that kind of material behavior must 
involve the consideration of the cooperative behavior or 
the self-organization of disordered structures. 

In this study we will try another route to tackle disor- 
dered structures. Before further discussing our approach, 
we need to take a close look at how disordered structures 
are formed under thermal fluctuations. It is believed that 
a crystalline material has a purely ordered structure at 
absolute zero and a completely disordered structure near 
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(a) (b) (c) 

FIG. 1. Schematic representation of the emergence of dis- 
torted crystal lattice; (a) the normal crystal lattice without 
distortion; (b) the distorted crystal lattice; (c) another ver- 
sion of the distorted crystal lattice. Obviously, both (b) and 
(c) have an equal probability of occurring. 



its melting point; at temperatures far below the melt- 
ing point, the considered material can by no means be 
assumed to be perfectly ordered. As diagrammatically 
shown in Fig. [1], there always exists the probability 
that certain atoms in the normal crystal lattice, which 
is shown in Fig. [la], could gain extra kinetic energy to 
move quasi-permanently away from their original equi- 
librium positions and distort crystal lattice as shown in 
Figs, [lb] and [lc] due to thermal fluctuations. The lo- 
cal disordered structures corresponding to such distorted 
crystal lattice arc defined as high-temperature structure 
precursors (HTSPs). At new equilibrium positions, the 
distorted lattice possesses higher potential energy due to 
the induced local strain. Thus, HTSPs associated with 
such lattice distortion are often metastable: they could 
disappear or even "/iop" to other locations under thermal 
fluctuations. Since the structures shown in Figs, [lb] and 
[lc] have an equal probability of occurring, they can also 
switch from one to the other in either direction under 
external perturbations. 

We now consider whether HTSPs could behave cooper- 
atively. Since the formation of HTSPs is caused by ther- 
mal fluctuations, their quantity is proportional to tem- 
perature. Therefore, as temperature rises, the quantity 
also increases. At a certain temperature Tjv, the quantity 
of HTSPs reaches a threshold so that they will start to 
interact with each other, which can be regarded as their 
cooperative behavior or self-organization, and then form 
a unique ncmatic phase. The driving force behind this 
structural transformation is the competition between en- 
ergy and entropy. Let's consider a single-phase metallic 
material containing HTSPs. Its Gibbs free energy in the 
absence of external fields can be written as G = U — TS; 
clearly, the internal energy U is increased due to the in- 
crement of the local strain energy generated by the dis- 
torted lattice shown in Fig. [lb] or [lc]. If temperature T 
remains unchanged or changes slowly, the entropy S must 
rise to reduce G. The simplest way to increase S is that 
the chemical bonds between the atoms at the distorted 
lattice points are partially broken so that the correspond- 
ing disordered structures could gain more freedom to ro- 
tate and are oriented along local preferred directions to 
form a nematic phase, which leads to a decrease in the 
orientational entropy but an increase in the positional en- 
tropy and, eventually, results in a net increase in the total 



entropy. This kind of structural transformation is very 
common in the first order nematic-isotropic phase tran- 
sition in liquid crystals [H, [H| . For the distorted lattice 
shown in Fig. [lb] or [lc], its distortion direction has an 
equal probability of pointing parallel or anti-parallel to 
a given direction (for instance, one of vertical directions 
shown in Fig. [lb] or [lc]). If wc use v to represent this 
given direction and approximately treat those HTSPs as 
molecules, then their orientational order parameter, S op , 
can be written as [l4[ 

Sop = 1(3 (v\n) 2 - 1) = i((3cos 2 ^ - l)), (1) 

where ( ) represents the average; v l is defined as the given 
direction of the disordered structure located at the posi- 
tion i; ft is usually called the director that represents a 
particular direction; 6 l is defined as the angle between v l 
and ft at the position i. If < S op < 1, we can say that 
HTSPs cooperatively form a nematic phase. S op = 1 
corresponds to an ideal case, in which all HTSPs are 
perfectly aligned. If temperature continues to rise, at a 
certain point T = TV/, the corresponding thermal en- 
ergy will be large enough to disturb HTSPs, which forces 
HTSPs to be randomly oriented and S op = 0. There- 
fore, when T > Tni, all HTSPs behave like a normal 
liquid, which has an isotropic phase. In this paper, we 
only consider the case, in which Tjv < T < Tjvj. 

In the above discussion, we use an orientational order 
parameter to represent, in the statistical sense, the co- 
operative behavior of HTSPs. Before giving further dis- 
cussion, we have to briefly introduce the concept of the 
order parameter. In 1937, Landau introduced a specific 
thermodynamic variable, i.e., the order parameter, to de- 
scribe his idea of broken symmetry in continuous phase 
transitions [l6|]. In the absence of the external field and 
in the vicinity of the critical point, T c , the Landau free 
energy F of a unit-volume material can be written as 
follows [3 

F = F + \p 2 + h)P A + • • • , (2) 

where P is defined as the order parameter and Fq is 
the free energy that is independent of P; both a and 
b are coefficients; Within the framework of the Landau 
theory, P has the following values: P = if T > T c , 
which corresponds to one phase with higher symmetry, 
and P = ± constant if T < T c , which corresponds to an- 
other phase with lower symmetry. If temperature T con- 
tinually changes, we can say that the material undergoes 
a continuous phase transition when T passes through the 
critical point. In this study, we generalize Landau's con- 
cept of the order parameter to represent the collective 
response of a material when disturbed by external fields. 
If a uniform stress a is applied to a unit- volume metal- 
lic material, the collective response, or strain e, of the 
material can be defined as its order parameter. We now 
consider the behavior of the above-defined nematic phase 
under the perturbation of a. Since HTSPs are usually 



3 



metastable, a, just like thermal fluctuations, will disturb 
them. According to Le Chatelier's principle, the nematic 
phase as the ensemble of HTSPs would then undergo a 
specific structural change, in which HTSPs would tend to 
huddle together and grow up, to counteract any imposed 
deformation by a |l5|. Therefore, the effective macro- 
scopic stress inside the material will be a e f f = a — aedis 
here the negative sign is due to Le Chatelier's princi- 
ple, tdis is the effective strain induced by HTSPs in the 
nematic phase, and a is the coefficient related to the 
nematic phase. For convenience, we can further write 
<J e ff = <J — aedis = (1 — k)a and define k as 



k 



koS op (T - T i 



N) 



c iV 



(T N <T< T NI ), 



(3) 



where both k and kg are dimensionless coefficients. Then 
the Gibbs free energy of that unit-volume metallic mate- 
rial can be written as 

G = U - TS - (1 - k)ae = U-TS-(1- k)Ee 2 , (4) 

where one can see that k actually represents, in the statis- 
tical sense, the fraction of the potential energy, which is 
generated by the cooperative movement of HTSPs. Now 
we can briefly summarize our approach: both the crys- 
talline phase and the nematic phase as the ensemble of 
HTSPs in a test specimen are integrated together by us- 
ing a modified mean field theory (MMFT), in which the 
former corresponds to the collective response of the crys- 
talline phase to external stimuli and the latter the ex- 
tremely slow counteractive and fluctuating response of 
the nematic phase of HTSPs in the specimen. 

If we assume that a is a function of time, then, as 
time goes, k will increase due to Le Chatelier's principle. 
Thus the modulated strain yl — ke is also a time- varying 
value. Clearly, for a perfect single-phase crystalline ma- 
terial, k w at the beginning when the external force 
is applied. At this moment, it is safe to write down the 
relationship a — Ee if the value of a is small. How- 
ever, k will continue to increase under the perturbation 
of the time-varying external force. In order to deter- 
mine how the evolution of k could eventually alter the 
linear stress-strain relationship, we need to exploit the 
Landau-Khalatnikov equation to investigate the relax- 
ation behavior of the modulated strain. 



III. FRACTIONAL POWER-LAW MATERIAL 
BEHAVIOR IN CHARACTERISTIC 
RELAXATIONS 

In 1954, Landau and Khalatnikov formulated an im- 
portant equation, which is written below fl7j . 
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dP 
dt 



dF 

ap' 



(5) 



where 7 is a kinetic coefficient and is considered to be in- 
dependent of temperature (l3l . [l8j ; t represents the time 



variable; the definitions of both F and P have been given 
in Eq. @. This equation is usually called the Landau- 
Khalatnikov equation that was originally developed to 
describe the critical slowing down of the fluctuation of the 
order parameter on approaching the critical point [l8j]. 
But it has a more profound physical significance; it can be 
interpreted as another version of the law of conservation 
of energy, i.e., the kinetic energy associated with the fluc- 
tuation of the order parameter and dissipated during the 
corresponding relaxation process is equal to the decrease 
in the Landau free energy. Therefore, this equation can 
be exploited to analyze the relaxation behavior of the or- 
der parameter in general cases. In this study, we replace 
Landau's concept of the order parameter (the sponta- 
neous strain) by a generalized order parameter (the in- 
duced strain). The uniqueness of our approach is that 
both the above-mentioned collective and slow counterac- 
tive responses can be integrated together in the mathe- 
matical expression of the Gibbs free energy to represent 
the competition between the crystalline phase and the 
nematic phase in the test specimen. Therefore, we can 
investigate this competition or, in other words, how the 
time-varying k could eventually alter the linear stress- 
strain relationship via relaxation processes. 

We now need to modify the Landau-Khalatnikov equa- 
tion for our studies. Notice that \aP 2 does positive con- 
tribution to the Landau free energy given in Eq. ^ 
but (1 — k)Ee 2 does negative contribution to the Gibbs 
free energy given in Eq. (j4]), so we have to modify the 
Landau-Khalatnikov equation as follows. 
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de 
df 



dG 



(6) 



Substitute G given in Eq. (jj) into the above equation, 
we get the following result. 



de dG . 



(7) 



From this equation, we can see that the SI unit of 7 is 
N • s/m 2 , which is similar to that of dynamic viscosity (in 
fact, it has been pointed out that 7 is viscosity for liquid 
crystals [HI). Solve this equation, we have, 



e s exp 



2(1 - k)E 

7 



(8) 



where e s = e|t=o, which can be interpreted as the induced 
strain under static deformation. For simplicity, we call 
e s the static strain. Eq. ([5]) represents the relaxation 
behavior of the induced strain after an instantaneous ex- 
ternal perturbation. Since the relaxation process consid- 
ered here is ultimately connected to the motion of atoms 
in crystalline lattices or partially ordered structures, we 
would introduce an atomic relaxation time, r, to charac- 
terize the relaxation behavior at the atomic level. Now 
we can approximately estimate the magnitude of the pa- 
rameters given in Eq. ((5J. For a fully relaxed process, we 
could regard t as t sa r; since r ~ ^L, where fn is defined 
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as the Debye frequency that is about fn ~ 10 13 Hz [l9j |. 
so t ~ lCP 13 s; for most metallic materials, E ranges from 
10 9 N/m 2 to 10 n N/m 2 ; k ranges from to 0.5 (we will 
discuss this later); since 7 represents the kinetic behav- 
ior at the atomic level, it is reasonable to assume that 
7> L Actually, 7 could be regarded as the viscosity, 
which describes the viscous behavior of the ncmatic phase 
of HTSPs; therefore, it must be very large. Thus we can 
estimate that | - 2(1 ~ fe)E £| < 1. By using Eq. flAT) in 
Appendix, we linearize Eq. (j8|) as follows 

2(k- l)E ' 
e«e s l + >—t . (9) 

7 

The above equation only represents a single relaxation 
after an instantaneous external perturbation; if we want 
to use it to investigate fatigue phenomena, the special 
cases of the fractional power-law material behavior, we 
must find a series of such atomic relaxations to form a 
continuous and progressive relaxation process in the test 
specimen when it is subject to cyclic loading. For this 
purpose, we consider three loading signals that are shown 
schematically in Fig. [2] . The first signal S-I, one period 
or cycle of it is shown in Fig. [2a], is a one-direction 
loading signal; there is no reversal in such a signal so 
that its cycle number is usually used to demonstrate fa- 
tigue phenomena. If we define T p as the period of S-I, 
in the first half of T p , a constant load is applied; no load 
in the second half. If T p r, the atomic relaxation will 
fully occur once at -f. The second loading signal S-II, 
one period or cycle of it is shown in Fig. [2b], is a fully 
reversed loading signal; there are two reversals in one cy- 
cle of S-II and its number of reversals is usually used to 
demonstrate fatigue phenomena. If we also define T p as 
the period of S-II, in the first half of T p , a constant load 
is applied; in the second half, the same load but in the 
opposite direction is applied. If T p 3> t, the atomic re- 
laxation will fully occur twice at both and T p , respec- 
tively. The third loading signal S-III, shown in Fig. [2c], 
is a constant loading signal; there is no reversal in S-III 
and the atomic relaxation does not occur during the time 
in which S-III is applied. For simplicity, we will mainly 
choose S-I with N cycles as the applied loading signal in 
our studies. Since both S-I and S-II are commonly used 
in fatigue studies in practice, we will later convert our 
derived formulas to the ones based on S-II via e — > e/2 
and N — > 2N. In this way, we will be able to present 
the fatigue relationships connected to both S-I and S-II, 
respectively. Here we have to emphasize that the fatigue 
phenomena connected to S-I might be different from the 
ones connected to S-II. The reason for this can be ex- 
plained in such a way that the value of k increases more 
quickly under the S-II loading test than that under the 
S-I loading test. 

Because the applied loading signal S-I has N cycles, t 
given in Eq. (|9|) should be re-counted as t = Nr. Substi- 
tute it into Eq. ([S]), we have the following relationship. 

e « e s [1 + (k - 1) P tN] , (10) 



where p = — . This equation represents the aforemen- 
tioned continuous and progressive relaxation process oc- 
curring in the test specimen when it undergoes a cyclic 
loading test. We now try to use it to derive different 
fatigue relationships. If assuming \prN\ <C 1, then, ac- 
cording to Eq. (|B1[) in Appendix, we can further simplify 
Eq. (fT0|) and get a new result as follows 

e^e s {l+(yrN) k - 1 =e s (^+pr\ N^ 1 
= e sm N-^ = e sm N-P, (11) 

here e sm = e 8 (ir + pr) k 1 . Making the following con- 
version, e — > e/2 and N — > 2N, we get the same relation- 
ship but based on S-II, which is given below. 

I « e sm (2N)-^ = e sm (2N)-P. (12) 

This equation describes a general relationship between e 
and N; here N may not be the fatigue life Nf. When 
N < Nf, both e sm and (3 are functions of N. Notice that 
we did not mention the amplitudes of the loading signals 
defined in Fig. [2]. Therefore, e and its corresponding 
applied stress are chosen arbitrarily. Now, if we choose an 
applied stress that is large enough for plastic deformation 
to occur in the test specimen, a fixed fatigue life Nf can 
be found to correspond to the induced plastic strain e p . 
Then we can get the following equation 

^Kef(2Nf)-(i-V=e f (2Nf)-e, (13) 

where e/ = e sm \N=Nf, which is defined as the ductility 
coefficient. This equation is the mathematical expression 
of the Coffin-Manson law, which was proposed indepen- 
dently by L. F. Coffin in 1954 and S. S. Manson in 1953 
Q to describe an empirical relationship between the in- 
duced plastic deformation and the number of reversals to 
failure. j3 = 1 — k given in the above equation is defined 
as the Coffin-Manson exponent. It has been observed 
that this exponent possesses a remarkable universality; 
P ~ 0.5 has been found in single-phased metallic materi- 
als whatever crystalline or polycrystallinc structures they 
have [1, [2(| • It is easy to see that the Coffin-Manson law 
is only the special case (when N = Nf and k ~ 0.5) of 

Eq. eg. 

Unfortunately, at this moment, there is no convincing 
theory that can be used to explain the above-mentioned 
universal law observed in the Coffin-Manson exponent 
Q . We now try to interpret the physical meaning of this 
exponent within the framework of the MMFT. As we 
have already discussed, the effective macroscopic stress 
inside the test specimen is (Jeff = o — atdis = (1 — k)a. 
It is well known that there is an essential prerequisite 
for the rationality and validity of any thermodynamic 
model, i.e, the specimen should be a continuum at any 
time when it is subject to an external perturbation. In 
thermodynamics, material properties are studied by us- 
ing continuous functions so that the considered material 
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must be treated as a continuum [l4|. Now we need to 
take a close look at the meaning of /3 = 1 — k ~ 0.5 or 
k ~ 0.5. Under the framework of the MMFT, for a single 
phase crystalline material, the value of k characterizes 
the counteractive response of the nematic phase as the 
ensemble of HTSPs to external perturbations; whereas 
that of 1 — k represents the collective response of the 
crystalline phase to external perturbations. Since the 
chemical bonds between certain atoms of HTSPs are as- 
sumed to be partially broken, they cannot be considered 
the integral part of the crystalline lattice of the test spec- 
imen anymore. Therefore, k represents, in the statistical 
sense, the fraction of the total atoms that are not the in- 
tegral part of the crystalline lattice in the test specimen. 
In other words, k = 0.5 represents a critical point; below 
this point (k < 0.5 or f3 > 0.5), the specimen is still a 
continuum; above this point (k > 0.5 or f3 < 0.5), the 
specimen is not a continuum anymore. Thus, from the 
viewpoint of thermodynamics, the Cofhn-Manson law is 
only valid when k < 0.5. It might also be interesting 
to interpret the meaning of /3 ~ 0.5 from the viewpoint 
of deformation and fracture of materials. Under a small 
external perturbation, the crystalline lattice will undergo 
harmonic vibration, which corresponds to the elastic de- 
formation, and HTSPs will "flow" , which corresponds to 
the plastic deformation, in the test specimen. If the per- 
turbation is large enough, both the crystalline lattice and 
HTSPs will undergo plastic deformation simultaneously 
in the test specimen. When k < 0.5 or f3 > 0.5, the 
induced clastic deformation dominates; thus we can say 
that the global deformation in the test specimen, in the 
statistical sense, is elastic. Of course, the local plastic de- 
formation due to the movement of HTSPs could also be 
induced in the test specimen even when k is very small; 
but this cannot alter the fact that the global elastic de- 
formation dominates under such circumstances. When 
k > 0.5 or j3 < 0.5, however, the induced plastic de- 
formation dominates; thus we can say that the global 
deformation in the test specimen, in the statistical sense, 
is plastic. Therefore, when (3 = 1 — k —> 0.5, the global 
plastic deformation will occur and the fracture crack will 
emerge and start to grow in the test specimen. In this 
sense, fatigue can be regarded as the precursor to frac- 
ture. This also explains why the global plastic defor- 
mation can be observed in the Coffin-Manson law when 
= 1 — k — > 0.5 in single phase crystalline materials. In 
practice, since there always exist defects in real metal- 
lic materials and/or alloys and these built-in disordered 
structures could lower the value of k, j3 usually ranges 
from 0.5 to 0.7 @]. Clearly, the Coffin-Manson law is one 
of the fractional power-law relationships, in which k ^ 
and f3 = 1 — k < 1. 

The fatigue phenomena described by the Coffin- 
Manson law are usually defined as low-cycle fatigue 
(LCF), in which the amplitude of loading signals is large 
and the number of cycles to failure is small. There is an- 
other type of fatigue, high-cycle fatigue (HCF), in which 
the amplitude of loading signals is very small and the 



number of cycles to failure is very large. We now try 
to derive one of important HCF laws, the Basquin equa- 
tion, which was proposed by O. H. Basquin in 1910 [j| 
to describe an empirical relationship between the elas- 
tic strain amplitude and a large number of reversals to 
failure. Since the number of cycles to failure in HCF is 
usually more than 10 4 0, we define the cycle number 
used in the following derivation as Nh so that it can be 
distinguished from N used previously. Similarly, we let 
t = Nut and then substitute it into Eq. ©. We have 
the following equation 



e s [1 + (k - l)prN h 
l + r](k 



(14) 



here rj is a coefficient, which is less than 1; we will explain 
why we need this coefficient later. The above equation is 
the corresponding continuous and progressive relaxation 
process occurring in the test specimen when it undergoes 
a HCF test. If assuming \prNh\ <C 1, then, according to 
Eq. (|B1I) in Appendix, we can further simplify the above 
equation and get a new formula as follows 



Nh 

1 + PT — 

V 



T,(k-1) 



1 

N~h 



,N 
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r,(k-l) 



Q~sn - v (l-k) 

E h 



E h ' 



V(k-l) 



where c. 



(A. 

\N h 



(15) 

k) is an ex- 



^N h ' T) ) ' 

poncnt. Since the deformation involved in HCF is of- 
ten believed to be mainly elastic, e sn given in the above 
equation is usually replaced by ^g 2 -. Making the follow- 
ing conversion, e — > e/2 and Nh — > 2N h , we get the same 
formula but based on S-II, which is written as follows 



E 



■{2N h y 



- V (l-k) 



E 



<2N h y 



(16) 



This equation describes a general relationship between 
e and N^, similarly, Nh may not be the fatigue life Nf 
and both a sn and b are functions of Nh when Nh < Nf. 
Now, if we choose an applied stress that can only induce 
clastic deformation in the test specimen at the initial 
stage, a long fatigue life Nf can be found to correspond 
to the induced strain e e . Then we can write the following 
equation 



1L 

E 



-V(l-k) = Zl 

E 



(17) 



where Of = cr sn |jv h= jv / , which is defined as the fatigue 
strength coefficient. This equation is the mathematical 
expression of the Basquin equation and b = rj(l — k) 
is defined as the Basquin exponent for convenience. In 
practice, b ranges from 0.125 to 0.2 Q. It is clear that the 
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Basquin equation is only the special case (when Nh = Nf 
and k ~ 0.5) of Eq. p^|) . It is also clear that the Basquin 
equation is one of the fractional power-law relationships. 



J Tp 






-0 

Tp/2 












(b) 



(c) 



FIG. 2. Schematic representation of loading signals; (a) S-I: 
the characteristic relaxation occurs only once within one T p ; 
(b) S-II: the relaxation occurs twice within one T p ; (c) S-III: 
no relaxation occurs; here T p is the period of loading signals. 



From what have been derived above, it is easy to 
see that there is no fundamental difference between the 
Coffin-Manson law and the Basquin equation. Both of 
them arise from the competition between the crystalline 
phase and the nematic phase as the ensemble of HTSPs 
during the aforementioned continuous and progressive re- 
laxation process in the test specimen; when approaching 
the numbers of reversals to failure in both cases, k is 
going to reach its limit, i.e., k — > 0.5. Only k varies dif- 
ferently in both cases; it increases quickly in the case of 
the Coffin-Manson law but slowly in that of the Basquin 
equation. Such dynamic behavior of k might cause prob- 
lems in the latter case. If we assume that b = 1 — k; since 
k may vary slowly in a HCF test, there always exists the 
possibility that when b = 1 — k is still around 1 (k ~ 0) 
whereas Nh has become very large but has not reached 
Nf, which will result in 4f —> 0. In order to prevent such 
a problem from occurring, we have to introduce a coeffi- 
cient rj into the Basquin equation to reduce the value of 
b. Obviously, for b = rj(l — k), rj must be less than 1. But 
the choice of the value of rj is quite arbitrary, which can 
also cause problems for HCF studies in practice. Proba- 
bly, this is one of the reasons that the Basquin equation, 
even widely used, is not so accurate as the Wohler equa- 
tion or the Wohler curve Q , another important HCF law 
proposed by A. Wohler in 1870 @,0|- 

We now try to derive the Wohler equation. Assum- 
ing that \(k - l)prN h \ < 1, then, using Eq. (|A2p in 
Appendix, we can modify Eq. (|14[) and then get a new 
equation as follows 



e«6,-E>[l- (k-l)pTN h ] 
1 



1 - In 

logio(JVfc) 



tt + (1 - k)pr 



logio(e) 
e s tp - 2.3e s logi (iV/ l ), 



(18) 



where ip = |l — In + (1 — k)pr |. Since the above 
equation represents a HCF relationship, the deformation 



involved in this equation is usually defined as the elas- 
tic one. Thus we can multiply both sides of the above 
equation by E and get the following result 



otw - P w \og w {N h ), 



(19) 



where a w = eE, a w = e s iJjE, and j3 w = 2.3e s E. Simi- 
larly, this equation describes a general logarithmic rela- 
tionship between a w and Nh', here Nh may not be the 
fatigue life Nf. Now, if we choose an applied stress a w 
that can only induce elastic deformation in the test spec- 
imen at the initial stage, a long fatigue life Nf can be 
found to correspond to a w . Then we can get the follow- 
ing result 



tf - p w log 10 (N f ), 



(20) 



where oef = a w \N h =N f ■ This equation is the mathemat- 
ical expression of the Wohler equation [(| 0]. One can 
see that the Wohler equation is actually the special case 
(when N h = N f and k ~ 0.5) of Eq. (JTj5J). 

In the Basquin and the Wohler equations, both k ~ 0.5 
and Nh = Nf (for perfect single-phase metallic materials) 
demonstrate that the co-existence of the induced elastic 
deformation and the induced plastic deformation in the 
test specimen has reached the balance point, which con- 
tradicts the widely used approximation method that the 
induced deformation in both cases is treated as the elas- 
tic one so that the stress variables instead of the strain 
variables are used in both equations. In practice, such a 
treatment is particularly straightforward and convenient; 
but, in theory, it may cause some problems since the in- 
duced deformation is by no means purely elastic under 
those circumstances. 

It might also be worth considering a special case, in 
which the loading signal is S-III shown in Fig. [2c] . There 
is no reversal in S-III and the atomic relaxation does 
not occur during the time in which S-III is applied. If 
we assume that the amplitude of S-III is large enough 
to induce plastic deformation in the test specimen, the 
value of k could also reach 0.5 under static or quasi-static 
deformation. Since the relaxation does not occur under 
such deformation, without loss of generality, we let t = r 
and substitute it into Eq. ([§]); the result is written as 
follows. 



e s exp 



-2(1 -k)E- 
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e,exp[-(l-k)pr], (21) 



here pr is dimensionless; e represents the induced final 
plastic strain and e s can be regarded as the initial static 
strain. As we have already discussed, we know that | — 
(1 — k)pr\ <C 1. By taking advantage of Eq. (| Al|) in 
Appendix, we can simplify the above equation and get 
the following formula 



or 



(1 - k)pre s 



e + (1 - k)pre s 



(22) 



(23) 
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Since pr < 1, by using Eq. (|B1I) in Appendix, we can 
further simplify the above equation and get a new result 
as follows 



1 



PT 



-1-fe 



(24) 



This equation is the fundamental equation that can be 
interpreted as follows. During the time, in which S-III is 
applied, the value of k increases as time goes. The above 
equation demonstrates the relationship between the ini- 
tial static strain e s and the final plastic strain e; for a 
perfect single phase crystalline material, k = at the 
initial stage and k —> 0.5 at the final stage. Now we try 
to alter the mathematical expression of the above equa- 
tion. Simply multiplying both sides of the above equation 
by an effective modulus, E, we get 



K sh e r a 



(25) 



where a = (e s + l)E, a y = eE, and n = 1 — k; 

f \ 1-fe _ 

K s h — I j- + pr J E. This equation is the mathemat- 
ical expression of the Ludwik's equation, which was pro- 



posed by P. Ludwik in 1909 [2lJ. Here we need to give 
extra explanation of this equation. In practice, the time 
interval between the initial stage (k = 0) and the final 
stage (k — > 0.5) is short since S-III is large, therefore, 
we have the following approximate definitions: a is the 
stress, <j y represents the yield stress, and e s is regarded 
as the plastic strain; n — 1 — fc is usually called the strain 
hardening exponent. If we re-write Eq. ([24]) as follows, 



1 = e 



1 

h pr 



i-fe 



,i-fc 



(26) 



here e approximately represents the induced plastic strain 
at any given time between the initial stage and the fi- 
nal stage. Similarly, multiplying both sides of the above 
equation by E, we have 



K sh e 



l-fc 



K sh e r s 



(27) 



here a = eE. This equation is the mathematical expres- 
sion of the Hollomon's equation, which was proposed by 
J. H. Hollomon in 1945 In practice, for most metallic 
materials and/or alloys, n generally varies between 0.2 
and 0.5 Q. The fractional power-law relationships be- 
tween the applied stress and the induced strain described 
by both the Ludwik's and the Hollomon's equations are 
usually defined as the strain hardening effect. The clas- 
sical explanation of this effect is that the induced plastic 
deformation would cause the generation of additional dis- 
locations in the test specimen; the more dislocations, the 
more they would likely become pinned together, which 
reduces the mobility of dislocations and prevents further 
deformation from occurring ||. Obviously, this effect 
can also be explained by the model proposed here, i.e., 
the material hardening or strengthening is due to the 



fact that the aforementioned disordered structures would 
tend to huddle together and grow up, due to Le Chate- 
liers principle, to counteract any imposed deformation by 
the applied stress. 

It is interesting to note that the strain hardening expo- 
nent n ranges from 0.2 to 0.5, which corresponds to the 
value of k varying from 0.5 to 0.8. As we have discussed, 
k < 0.5 is the prerequisite for the validity of any ther- 
modynamic model. Why is k greater than 0.5 in both 
the Ludwik's and the Hollomon's equations? We try to 
explain this as follows. In order to derive these two equa- 
tions, we multiply both sides of Eq. (f2~4")l by E. In doing 
so, we actually assume that there are linear relationships 
between a and (e s + k), between a y and e, etc. Ob- 
viously, this simple approach cannot result in perfectly 
correct mathematical expressions since the deformation 
involved here is plastic. This might be one of the reasons 
that the measured n could be less than 0.5 in practice. 
Theoretically, compared with both the Ludwik's and the 
Hollomon's equations, Eq. (|24p may be a better model 
for the mathematical description of the strain hardening 
effect. 



IV. CONCLUDING REMARKS 



In this study, the origin of the fractional power-law 
material behavior has been investigated. We generalize 
Landau's concept of the order parameter to represent the 
collective response of the crystalline phase of a test spec- 
imen when disturbed by external fields and treat its HT- 
SPs as a nematic phase or a partially ordered liquid over 
a wide range of temperature. Both responses of the crys- 
talline and the nematic phases to external fields are then 
integrated together within the framework of the MMFT. 
In doing so, the relaxation behavior of the generalized 
order parameter of the specimen can be studied via a 
modified Landau-Khalatnikov equation. In view of what 
have been derived and discussed above, one can see that 
the fractional power-law material behavior is closely re- 
lated to the nematic phase of HTSPs; furthermore, we 
can conclude that, as the special cases of the fractional 
power-law material behavior, different fatigue phenom- 
ena are actually governed by the competition between 
the crystalline phase and the nematic phase as the en- 
semble of HTSPs during a series of atomic relaxations in 
the test specimen. 
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Appendix A: Series expansions of exponential and 
natural logarithm functions 



Similarly, if \x\ <C 1, the above equation can be simpli- 
fied as follows, 



For an exponential function e x , its Maclaurin series 
expansion can be written as 



In (1 - x) 



(A2) 



n=0 



X" 



= 1 + x + _ + _ + ... . (for all x) 

If | a; | <C 1, the above equation can be simplified as 
follows, 



e x w 1 + x. 



(Al) 



For a natural logarithm function In y = ln(l — x), its 
Maclaurin series expansion can be written as 



n=l 
= —X — ■ 



X" 

n 



(-1 < x < 1) 



Appendix B: Binomial series expansion 

For a function f(x) = (1 + its binomial series 
expansion can be written as 

n\ 

For simplicity, we here only consider the case that £ 
is a real number. If |x| < 1, this expansion converges 
absolutely for any number £. If \x\ <C 1, we can get the 
following approximation formula, 



{l + x) e m l + fcc. 



(Bl) 
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